arXiv:1507.02138vl [astro-ph.HE] 8 Jul 2015 


Mon. Not. R. Astron. Soc. 000,^?? (2015) Printed 9 July 2015 (MN style file v2.2) 

The magnetic field strengths of accreting millisecond 
pulsars 

Dipanjan Mukherjee^*, Peter Bult^, Michiel van der Klis^ and Dipankar Bhattacharya^ 

^Research School of Astronomy and Astrophysics, The Australian National University, Canberra, ACT 2611, Australia 
^ Anton Pannekoek Institute, University of Amsterdam, Science Park 90f, 1098 XH Amsterdam, The Netherlands 
^Inter University Centre for Astronomy and Astrophysics, Post Bag f, Pune fllOOl, India. 


9 July 2015 


ABSTRACT 

In this work we have estimated upper and lower limits to the strength of the magnetic 
dipole moment of all 14 accreting millisecond X-ray pulsars observed with the Rossi 
X-ray Timing Explorer {RXTE). For each source we searched the archival RXTE data 
for the highest and lowest flux levels with a significant detection of pulsations. We 
assume these flux levels to correspond to the closest and farthest location of the inner 
edge of the accretion disk at which channelled accretion takes place. By estimating 
the accretion rate from the observed luminosity at these two flux levels, we place 
upper and lower limits on the magnetic dipole moment of the neutron star, using 
assumptions from standard magnetospheric accretion theory. Finally, we discuss how 
our field strength estimates can be further improved as more information on these 
pulsars is obtained. 
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1 INTRODUCTION 

Magnetic field strengths of neutron stars span a wide range; 
from > 10^'* G in magnetars via ~ 10^^ G in most radio pul¬ 
sars down to ~ 10® G in millisecond pulsars. The millisecond 
pulsars are thought to attain their fast spin and low mag¬ 
netic field strength through accretion in low-mass X-ray bi¬ 


naries (|Alpar et al.|1982 Radhakrishnan Sz Srinivasan|1982 


|Bailes| 1989p . The discovery of accreting millisecond X-ray 
pulsars (AMXPs, Wijnands fc van der Klis|1998 1 and tran¬ 
sitional millisecond pulsars (Papitto et al.| 2013bl has lent 


strong support to this picture. However, a full understand¬ 
ing of the evolutionary scenario will only emerge through 
the comparison of both spin and magnetic field distribu¬ 
tions of the LMXB neutron stars with those of transitional 
and millisecond radio pulsars. In this paper we undertake 
the task of estimating the magnetic field strengths of the 
AMXP population in a consistent approach. 

Magnetic field strengths of neutron stars are estimated 
in a variety of ways. For non-accreting radio pulsars the field 
strengths are inferred from the spin-down rate due to elec¬ 
tromagnetic torque (e.g. |Ostriker fc Gunn| 196^ [Beskin et al.| 
1984 Spitkovsky|2006 1. For anomalous X-ray pulsars, rough 


estimates of the field can also be obtained by modelling 
the non-thermal X-ray spectra with cyclotron and magnetic 
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Compton scattering processes in the magnetosphere (Giiver 
et al. p007|[2008t . 


For accreting systems, the available methods for esti¬ 
mating the field strength are less robust. X-ray pulsars with 
high surface field strengths (ft 10^^ G), may show resonant 
electron cyclotron lines in their X-ray spectra, which give 
an estimate of the local field strength in the line formation 
region | Caballero fc Wilms||2012 K However, the location of 
this region is unclear, resulting in uncertainties in the mea¬ 
sured dipole moment. For lower surface field strengths, as 
for AMXPs, cyclotron resonances move out of the X-ray 
band, leaving only indirect methods for estimating the field 
strength. 

In this paper we estimate the magnetic field strength of 
AMXPs using X-ray observations obtained with the Rossi 
X-ray Timing Explorer (RXTE). We assume that the ac¬ 
cretion disk is truncated at an inner radius that depends 
on the magnetic field strength and the accretion rate. Be¬ 
low the truncation radius the disk is disrupted and matter 
is forced to move along the magnetic field to the magnetic 
polar caps of the neutron star, creating localised hotspots 


responsible for X-ray pulsations (see, e.g., Poutanen fc Be- 
loborodov|2006 1. The detection of X-ray pulsations is taken 


to imply ongoing magnetically channelled accretion onto 
the neutron star, such that the highest and lowest flux lev¬ 
els with detected pulsations identify the range of luminosi¬ 
ties (and hence accretion rates) over which such accretion 
occurs. Assuming a disk-magnetic field interaction model, 
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these measurements can then be used to constrain the sur¬ 
face dipole held strength of the neutron star. 

To calculate the dipole moment from the set of hux lev¬ 
els, we consider the Ghosh & Lamb (1978 19791 model of 
disk-magnetic held interaction, applied in a manner akin to 
that adopted by Psaltis & Chakrabarty (19991. We assume 
magnetically channelled to mean that the disk truncation 
radius is outside the neutron star surface, but smaller than 


the co-rotation radius (Pringle & Rees 1972 Illarionov fc 
Sunyaev||1975 1, keeping in mind that these radii only set an 


approximate scale for the system. The accretion disk may 
not extend all the way down to the neutron star surface 
and channelled accretion may persist for a disk truncated 


outside the co-rotation radius (Spruit & Taam 


paport et al. 2004 Bult & van der Klis 20151. How this 


1993 Rap- 


choice of radii affects our magnetic held strength estimates 
is discussed in section [TS] 

We plan the paper as follows: in Sec. we review the 
theory of accretion-disk/magnetic-held interaction, the de¬ 
tails of our held strength estimation method and X-ray anal¬ 
ysis; in Sec. we describe the outburst history of the con¬ 
sidered sources and the results of our analysis; in Sec. we 
discuss the uncertainties in our method and how they affect 
our results; and hnally in Sec.we compare our results with 
previous estimates of the magnetic held strength. Technical 
details of the timing and spectral analysis are presented in 
Appendix and respectively. 


2 MAGNETIC FIELD STRENGTH 
ESTIMATION METHOD 


2.1 Theory of disk-stellar magnetic field 
interaction 


For accretion in a steady state, the inner truncation ra¬ 
dius depends on the balance between magnetic and material 
stresses. Equating the torque from the magnetic stresses and 
the angula r momentum hux ( [Ghosh fc Lamb||1979| [Rappa-j 
port et ^[20041^ one hnds that 


d(Mr^nK) 

dr 


= B^Bpr 


( 1 ) 


where M is the accretion rate, r the distance from the com¬ 
pact object, Q.K the Keplerian angular velocity at r, and Bp 
and B^ the respective poloidal and toroidal components of 
the magnetic held. For simplicity we calculate the torques at 
the truncation radius considering the spin axis to be aligned 
to the magnetic dipole, with the accretion disk being per¬ 
pendicular to both axes. The toroidal held component is 
produced due to shearing of the poloidal helds. Its strength 
is an uncertain parameter as it depends on various poorly 
understood processes like turbulent diffusion and magnetic 
reconnection ( |Wang|[i995| ). Another uncertain quantity is 
the radial extent of the accretion disk, Ar, over which matter 
couples to the stellar magnetic held and is channelled away 
from the disk. By expressing these uncertain quantities with 
the boundary layer parameter, 7 s = {B^/Bp){Ar/rt), the 


^ Viscous stresses at the truncation point have been ignored 


truncation radius rt can be related to the poloidal magnetic 
held and hence the dipole moment as 

Bp = ^-^i^{gmm^')^'\-^^\ ( 2 ) 

where G is the gravitational constant and M the neutron 
star mass. Assuming a dipolar magnetic held, Bp{r) = 1 

with /i the magnetic dipole moment, the truncation radius 
is given as 


n = iT ( 




1/7 


\GMM^ ) 


(3) 


The truncation radius is related to the classical Alfven ra¬ 
dius ta as rt = 7 ^'^^ 2 ^/^rA, where 


/ M \ / M ' 

"" Viol® gs-V V1-4M0, 


12/7 

I 

-1/7 


(4) 


Here Bs = nfR^ is the magnetic held strength at the equator 
and Rs the neutron star radius. 

At higher accretion rates the truncation radius will be 
closer to the neutron star surface. In this work we assume the 
disk to extend all the way down to the neutron star surface 
at the highest accretion rate. At the lowest accretion rate we 
assume the truncation radius to be at the co-rotation radius. 
Both radii are rough approximations as true behaviour of 
the accretion disk truncation radius depends on uncertain 
aspects like disk/held coupling and the local magnetic held 
topology (jKulkarni fc Romanova 2013 D’Angelo & Spruit 


2012 2010| Romanova et al. 20081. We discuss the effects 
and limitations of these assumptions in Sec. |4.3| 

By identifying the highest and the lowest accretion rates 
with ongoing magnetic channelling (conhrmed by the detec¬ 
tion of pulsations), we estimate the magnetic held as out¬ 
lined below. 

For all sources we adopt the canonical neutron star mass 
of M = 1.4 Mq and radius Rs = 10 km. The value of 7 b is 
highly uncertain and depends on where the accretion disk is 
truncated. To be conservative we take 7 b to vary between a 
wide range of 0.01 — 1 (jPsaltis &: Ghakrabarty|[l999|). 


(i) Lower limit on /i,: At the highest accretion rates, 
to observe pulsations, the magnetic held must be at least 
high enough to truncate the accretion disk at or above the 
neutron star surface. Thus by setting rt = Rs we obtain the 
lower limit on the dipole moment as 

= (5) 


We assume the mass accretion rate can be estimated from 
the bolometric luminosity as L = GMM/Rs, and estimate L 
from the observed luminosity in the X-ray band by applying 
a bolometric correction factor {L = Cboi-Lx). The typical 
reported values of the bolometric correction factor have a 


range of Cboi ~ 1 — 2 (jGilfanov et al.||1998 

Galloway et al. 

to 

0 

0 

to 

Gampana et al.|20031 

Mighari & 4’ender|2006| [Gasella 

et al. 

to 

0 

0 

00 

Galloway et al. 

2008 

1 . The mass accretion rate 
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then follows as 


M = 10^® g s"^ 
"" \lAMeJ 


tboiLx 


1.87 X 1036 erg s 


Rs 


10 km 



( 6 ) 


Expressing the X-ray luminosity in terms of the observed X- 
ray flux and the assumed distance {Lx = AtvcTF), we obtain 
the lower limit on the magnetic moment 


= 9.36 X lO"'* G cm-® 


1/2 


10“'^6 erg cm“3 3-1 
M 


1/2 


d 

10 kpc 


1.4 Mq 


-1/4 / \ 6/4 


10 km 


(7) 


where Emax is the highest observed X-ray flnx with pulsation 
and we adopted the boundary values of 75 and esoi that 
give the lowest magnetic moment. For a detailed discussion 
of these assumptions we again refer to sec. 

(ii) Upper limit on p: We assume magnetic channelling 
to be centrifugally inhibited if the accretion disk is outside 
the co-rotation radius 


/ GM Y''® 
V(27ri2a)2y 


( 8 ) 


where Vs is the spin frequency. We then obtain an upper 
limit on the dipole moment by setting the truncation radius 
at the co-rotation radius (rt = rP) for the lowest accretion 
rate with detected pulsations 

Mmax = (9) 


By again substituting the expression for the mass accretion 
rate we obtain the upper limit on the magnetic dipole mo¬ 
ment 



where Emin is the lowest observed X-ray flux with pulsation 
and we adopted the boundary values of 7 s and esoi that 
maximize the magnetic moment. 


2.2 Data analysis 

We analysed all AMXP outbursts observed with RXTE. 
The general structure of our analysis is as follows; for each 
RXTE observation we estimate the Crab normalized X-ray 
flux and then search for the presence of pulsations. We se¬ 
lect the highest and lowest flux observations with detected 
pulsations, and from spectral analysis measure the source 
flux. Using these flux measurements and the best distance 
estimate from literature, we obtain limits on the magnetic 
dipole moment from eq. [7|andp^ These limits are expressed 
in terms of the magnetic dipole field strength at the equator 
(Bp = /r/r), assuming a 10 km radius. The detailed proce¬ 
dure is outlined below. 


2.2.1 Timing analysis 


We initially estimate the 2-16 keV Crab normalised X-ray 
flux from the 16-second time-resolution Standard-2 data (see 
e.g. van Straaten et al.||2003 for details). To search for pul¬ 
sations we consider the high time resolution GoodXenon or 
(122 fis) Event mode data of the same observation, selecting 
only those events in energy channels 5-37 (~ 2 — 16 keV), 
which usually provides an optimal signal to noise ratio 
for the pulsations. The data were barycentred using the 
FTOOLS task faxbary, which also applies the RXTE fine 
clock corrections, thus allowing for timing analysis at an ab¬ 
solute precision of ~ 4/ts (Rots et al. 20041. We then take 
each ~ 3000 s continuous light curve (as set by the RXTE 
orbit), and fold it on the pulsar ephemeris (see appendix]^ 
to construct a pulse profile. For each profile we measure the 
amplitude at the fundamental frequency and that of its sec¬ 
ond harmonic (Hartman et al.||2008|. 


In standard procedures (see, e.g, Patruno & Watts 


20131, a pulsation is usually said to be significant if the mea¬ 


sured amplitude exceeds a detection threshold. This thresh¬ 
old is set as the amplitude for which there is only a small 
probability e that one the of observations in an outburst ex¬ 
ceeds it by chance. For observed amplitudes higher than this 
threshold we have a high confidence C = 1 —e that pulsations 
are detected (C = 99.7%). We can then consider the flux es¬ 
timates associated with the significant pulse detections and 
straightforwardly select the observations of highest and low¬ 
est flux. 


This approach is very conservative, as it sets a small 
joint false-alarm probability of detection for the entire out¬ 
burst, in spite of the fact that we can be certain that pul¬ 
sation are present in most observations. At the low flux end 
of the outburst, where the detection significance decreases 
with the count rate, this may cause us to miss pulsations. 


To overcome this issue we first reduce the number of 
trials by comparing the observed total count rate of an ob¬ 
servation with the X-ray background as estimated with the 
FTOOLS task pcabackest ( Jahoda et al.|2d06 |. We then set 
a minimal count-rate threshold and reject all observation 
for which the pulse amplitude cannot be detected above 
the noise level assuming the expected source contribution 
is 100% modulated. 


We then select all observations of an outburst that do 
not have significant pulsations according to the procedure 
described above. If the pulsar emission is indeed not present 
in these observations, then the distribution of measured am¬ 
plitudes and phases should correspond to the expected dis¬ 
tribution of random noise, i.e. the phases should be uniform 
and the squared amplitude should follow a x®-distribution 
for two degrees of freedom. We compare the distributions 
using a KS-test, again with a 99.7% confidence level. If the 
data is not randomly distributed we take out the highest flux 
observation that has a significant pulse detection at the sin¬ 
gle trial level and whose phase is consistent with the expec¬ 
tation from the timing model, and iterate until the sample is 
consistent with being random. The last removed observation 
is then taken as the lowest observed flux with pulsations. 
We note that in practice the sensitivity of this iterative ap¬ 
proach is limited by the small number of observations in the 
tail of an outburst and only rarely yields a lower flux pulse 
detection than through the initial procedure outlined above. 
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2.2.2 Spectral analysis 

For the spectral analysis of selected observations we used 
HE AS OFT version 6.12 and the calibration database 
(CALDB). Spectra were extracted from the Standard-2 data 
following the standard procedures outlined in the RXTE 
cookbool0 The background was again estimated using the 
FTOOLS task pcabackest. A dead-time correction was ap¬ 
plied to the spectra following the prescription in the RXTE 
cookbool0 The spectral fits were done in the 3-20 keV en¬ 
ergy range with XSPEC version 12.7.1. 

From the measured flux we calculate the 3-20 keV X- 
ray luminosity and convert to the bolometric luminosity by 
multiplying with the correction factor Cboi. The bolometric 
correction factor for a source depends on its spectral state, 
which in turn varies with accretion rate. The typical range of 
the bolometric correction factor is reported to be tboi ~ 1 —2 


iGilfanov et al.||1998 |Galloway et al.||2002| 

Gampana et al. 

2003lJMigliari & Fender|20061 IGasella et al. 

20081 IGalloway 

et al. 12008 1. To be conservative we use Cboi 

= 1 when calcu- 


lating /imin and eboi = 2 when calculating pmax 


For many of the AMXPs considered here, there is con¬ 
siderable contaminating background emission in the ob¬ 
served X-ray flux, for instance from Galactic ridge emission. 
To estimate the background contamination, we also mea¬ 
sured the bolometric luminosity for an observation in the 
tail of the outburst, where pulsations were not present and 
the light curve has asymptotically levelled off to a constant 
value. We assume that in such a state, accretion has ceased 
and the observed flux is purely due to background emission. 

The details of the spectral fit parameters for different 
sources are presented in Appendix. 


3 RESULTS 

In the following sections we present the magnetic field esti¬ 
mates we obtained from the timing and spectral analysis of 
all AMXPs observed with RXTE. For each source we briefly 
describe the outburst history, distance estimates and discuss 
specific details of our analysis. All results are summarized 
in Tables [T] and m 



55026 55028 55030 55032 55034 

Time (MJD) 


Figure 1. Upper panel: Light curve of the 2009 outburst of Swift 
J1756.9—2508 normalised to mCrab. Lower panel: pulse ampli¬ 
tude of the fundamental (squares) and the second harmonic (open 
circles). The observations with significant detection (with 99.7% 
confidence limit) of pulsation are marked with black, observations 
without a significant detection of pulsations are shown in grey. 



0 ^^^^_Lj_,_1_^_1_,_^_ 

52400 52410 52420 52430 52440 52450 
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Figure 2. Light curve of the 2002 outburst of XTE J0929—314. 
Legends are same as in Fig. 


3.1 Swift J1756.9-2508 


Swift J1756.9-2508 was first discovered with Swift in June 
2007 (Krimm et al. 2007a|b I, and 182 Hz pulsations were 
found with follow-up RXTE observations (Markwardt et al. 


2007). The source showed a second outburst in July 2009 
( Patruno et al.|2009b 2010a|. 

We find both the highest and lowest flux with pulsa¬ 
tions to occur during the 2009 outburst (Fig. [^, with a 
detection of pulsation for the outburst peak luminosity at 
MJD 55026.1 and lowest flux detection on MJD 55032.5, 
just before the light curve decays to the background level. 
The background contribution was measured from the last 
observation on MJD 55037.0. 


^ http://heasarc.nasa.gov/docs/xte/recipes/cook_book. 
html 

http://heasarc.gsfc.nasa.gov/docs/xte/r6cipes/pca_ 
deadtime.html 


The distance to the source is not known, but consider¬ 


ing its close proximity to the galactic centre (Krimm et al. 
|2007a| ), we assume a distance of 8 kpc. We then obtain a 
magnetic field range of 1.8 x 10^ G < B < 2.4 x 10® G. 


3.2 XTE J0929-314 


XTE J0929-314 was discovered in April 2002 with RXTE, 
and the 185 Hz pulsations were immediately detected 
I Remillard|2002 Remillard et al. 20021. The source has been 
detected in outburst only once, with the light curve shown 
in Fig. 

XTE J0929-314 is significantly away from the galactic 
plane (galactic coordinates: 260.1°, 14.2°) and shows a low 
neutral hydrogen column density (~ 7.6 x 10®® cm~®,|juett 


et al.|[2003). The background contribution for this source is 


therefore negligibly low, and indeed could not be measured 
as the non-pulsating observations at the end of the outburst, 
which have too few counts to constrain the spectrum. 
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Table 1. Flux range with pulsations of the analysed AMXPs in order of ascending spin frequency. 


Maximum flux Minimum flux Background flux 


No 

Name 

Spin 

(Hz) 

MJD 

Flux 

( erg cm~^ s“^) 

MJD 

Flux 

( erg cm“^ s~^) 

MJD 

Flux 

( erg cm~^ s 

1 

Swift J1756.9-2508 

182 

55026.1 

6.30 X 

10-10 

55032.5 

1.99 X 

10-10 

55037.0 

4.07 X 

10-11 

2 

XTE J0929-314 

185 

52403.5 

4.42 X 

10-10 

52442.3 

6.64 X 

10-11 

- 

- 


3 

XTE J 1807.4-294 

191 

52697.6 

8.19 X 

10-10 

52808.7 

3.51 X 

10-10 

52816.8 

7.25 X 

10-11 

4 

NGC 6440 X-2 

206 

55073.1 

2.62 X 

10-10 

55873.3 

3.36 X 

10-11 

55823.2 

1.34 X 

10-11 

5 

IGR J17511-3057 

245 

55088.8 

8.65 X 

10-10 

55124.0 

1.00 X 

10-10 

55118.2 

6.96 X 

10-11 

6 

XTE J1814-338 

314 

52814.3 

4.41 X 

10-10 

52844.1 

6.00 X 

10-11 

52852.9 

1.00 X 

10-12 

7 

HETE J1900.1-2455 

377 

53559.5 

1.15 X 

10“® 

53573.8 

3.84 X 

10-10 

- 

- 


8 

SAX J1808.4-3658 

401 

52563.2 

1.85 X 

10-9 

50936.8 

2.82 X 

10-11 

50935.1 

1.21 X 

10-11 

9 

IGR J17498-2921 

401 

55789.6 

1.13 X 

10“® 

55818.3 

4.44 X 

10-10 

55826.4 

4.23 X 

10-10 

10 

XTE J1751-305 

435 

52368.7 

1.50 X 

10“® 

52377.5 

3.97 X 

10-10 

52380.7 

6.51 X 

10-11 

11 

SAX J1748.9-2021 

442 

55222.5 

4.13 X 

10“® 

52198.3 

2.96 X 

10-® 

55254.6 

1.88 X 

10-11 

12 

Swift J1749.4-2807 

518 

55300.9 

5.24 X 

10-10 

55306.7 

2.67 X 

10-10 

55307.7 

2.41 X 

10-10 

13 

Aql X-1 

550 

50882.0 

8.74 X 

10“® 

50882.0 

8.74 X 

10-® 

50939.8 

1.34 X 

10-11 

14 

IGR J00291-I-5934 

599 

53342.3 

9.70 X 

10-10 

53352.9 

1.09 X 

10-10 

53359.5 

5.76 X 

10-11 


For each source we give the minumum and maximum flux with pulsations. We also give the corresponding 3-20 keV background flux 
used in calculating the magnetic dipole moment and the MJD of the analysed observations. 


Table 2. The magnetic field strength estimates (equatorial surface field) of all considered AMXPs. 


No 

Name 

Spin 

(Hz) 

Distance 

(kpc) 

Background 

corrected 

-^min 

(10® G) 

Bmax 

(10® G) 

-^K,max 

(10® G) 

Bl 

(10® G) 

ref. 

1 

Swift J1756.9-2508 

182 

8 

y 

0.18 

24.1 

7.2 

0.4-9 

a 

2 

XTE J0929-314 

185 

6 

n 

0.12 

11.5 

3.4 

- 


3 

XTE J1807.4-294 

191 

4.4 

n 

0.11 

18.6 

5.4 

- 


4 

NGG 6440 X-2 

206 

8.2 

y 

0.12 

7.6 

2.1 

— 


5 

IGR J17511-3057 

245 

7 

n 

0.19 

11.8 

3.1 

- 


6 

XTE J1814-338 

314 

8 

y 

0.16 

7.8 

1.8 

<5-10 

b 

7 

HETE J1900.1-2455 

377 

5 

n 

0.16 

10.0 

2.1 

- 


8 

SAX J 1808.4-3658 

401 

3.5 

n 

0.14 

1.77 

0.36 

0.7- 1.5 

c 

9 

IGR J17498-2921 

401 

8 

n 

0.20 

16.0 

3.2 

- 


10 

XTE J1751-305 

435 

7 

y 

0.25 

11.0 

2.1 

3.3 - 4.7 

d 

11 

SAX J1748.9-2021 

442 

8.2 

y 

0.49 

37.8 

7.2 

- 


12 

Swift J1749.4-2807 

518 

6.7 

n 

0.11 

7.7 

1.4 

- 


13 

Aql X-1 

550 

5 

y 

0.44 

30.7 

5.3 

< 9 

e 

14 

IGR J00291-I-5934 

599 

3 

n 

0.085 

1.9 

0.31 

1.5 - 2.0 

f 


(a)|Patruno et ^ (|2010a|l, (b)|Watts et al.|||2008| |; |Papitto et al.| j2007| l; |Haskell Si Patruno] l |2011[ |, (c) |Patruno et aT] | |2012j |, (d) [Riggl^ 
|et al.| l |2011a[ |, Di Salvo Burderi] < |2003a| l, (f) |Patruno| ( |201o{ " ' 

The values Bmin and -Bmax correspond to field strengths estimated using eq. and eq. |10[ The values Bk are upper limits to the field 
strength computed using eq. |12| following the modified expression of truncation radius as obtained by |Kulkarni Sz Romano^ < |2013l l (see 
Sec. |4.3| |. Bl are field strength measurements from literature (see references). The background correction column indicates if the 
background estimate (Table [3 was used when calculating the upper limit on the magnetic field estimate (see, e.g, Sec. |3.3||. 


There are no good estimates of the distance to the 
source, with the only constraint claiming d > 6 kpc based on 
estimates of average accretion rate Galloway et al. (20021. 
Using this distance we obtain a magnetic field range of 
1.2 X 10^ G < B < 1.2 X 10® G. 


3.3 XTE J1807.4-294 


XTE J1807.4-294 was discovered in February 2003 (Mark- 


wardt et al. 20031 and the 191 Hz pulsations were im¬ 


mediately found with the RXTE observations. The source 
has been in outburst only once. We find the highest and 


lowest flux with pulsations to occur on MJD 52697.6 and 
MJD 52808.7, respectively (Fig. [^. 


We measured the background contribution on 
MJD 52816.8, but note that this background level is 
similar to the lowest flux with pulsations, such that the 
uncertainty in the background estimation is larger than 
the apparent source contribution. To be conservative we 
calculate the upper limit to the magnetic moment without 
adjusting for the background. This implies that we take the 
low flux observation as an upper limit to the true lowest 
flux at which pulsations are present. If the presence of 
pulsations can be established at a lower flux level, the upper 
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XTEJ1807.4-294 



Figure 3. Light curve of the 2003 outburst of XTE J1807.4-294. 
Legends are same as in Fig. ^ 


IGR J17511-3057 



55085 55095 55105 55115 55125 

Time (MJD) 


Figure 5. Light curve of the 2009 outburst of IGR J17511-3057. 
Legends are same as in Fig. ^ 


NGC 6440 X-2 



Time (MJD) 


Figure 4. Light curve of NGC 6440 X-2 from July, 2009 to 
November, 2011. Legends are same as in Fig. [T] 


limit will decrease, thus tightening the allowed magnetic 
field range. 

There are no well defined estimates for the distance to 
this source. Some authors assume the source is near the 
galactic centre ( Falanga et al.|2005a Campana et al.|2003 l, 
and take the distance to be ~ 8 kpc. Others, however, 
estimate the distance at ~ 4.4 kpc ( [Riggio et al.|[2008[ ), 
by comparing the observed flux to the accretion rate in¬ 
ferred from the pulse timing analysis. Lacking a more ro¬ 
bust estimation of the distance, we adopt a distance of 

4.4 kpc. We then arrive at a magnetic field range estimate 
of 1.1 X 10^ G < B < 1.9 X 10® G. 


monitoring of this globular cluster, the background level is 
well established between MJD 55700-55850. We estimate 
the background flux on MJD 55823.2. 

The distance to the cluster NGC 6440 is well con¬ 
strained to d = 8.2 kpc (Valenti et al. 20071, and gives a 
magnetic field range of 1.2 x 10' G < B < 7.6 x 10® G. 


3.5 IGR J17511-3057 


IGR J17511-3057 was discovered with INTEGRAL in 


September 2009 (Baldovin et al. 20091, with subsequent 


RXTE observations discovering 245 Hz pulsations (Mark- 
|wardt et al.|200'^ . 

The outburst light curve shows a notable flare after 
MJD 55110 (Fig. which is attributed to simultaneous 
activity of XTE J1751-305 ( [Falanga et al.||201ip . Pulsa¬ 
tions are observed throughout the outburst, and reoccur on 
MJD 55124.0, some 10 days after the source appears to have 
reached the background level. We select this observation as 
the lowest flux observation with pulsations, and select the 
observation on MJD 55118.2 for the background flux. Since 
the flux difference between these observations is very small 
(Table we consider the low flux observation to be back¬ 
ground dominated and, like in Sec. |3.3| neglect the back¬ 
ground measurement in calculating the upper limit on the 
magnetic moment to obtain a more conservative estimate of 
the magnetic held strength. 

The distance to this source is estimated at <7 kpc, 
derived by assuming the type I X-ray bursts are Edding¬ 


ton limited (Altamirano et al. 2010b Papitto et al. 20101. 


Adopting this distance we obtain a magnetic field range of 
1.9 X 10^ G < B < 1.2 X 10® G. 


3.4 NGC 6440 X-2 

NGG 6440 X-2 is located in the globular cluster NGG 6440 
and was detected serendipitously with Chandra in July 2009 


(Heinke & Budac 2009 Heinke et al. 20101. Pulsations at 


206 Hz were discovered from subsequent RXTE observations 
( Altamirano et al.|2009 2010al. 

The outburst behaviour of NGC 6440 X-2 is atypical, as 
it shows brief outbursts of a few days with a recurrence time 
as short as one month (see Fig.|^. Due to the high-cadence 


3.6 XTE J1814-338 

XTE J1814-338 was discovered in June 2003 with RXTE 


and immediately recognized as a 314 Hz pulsar (Markwardt 
&: Swank|2003l. It has been detected in outburst only once. 


and shows pulsations throughout its outburst (Fig. [^. 

We measured the background flux from the observation 
on MJD 52852.9, one of the last RXTE observations of the 
outburst. 
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XTEJ1814-338 



Figure 6. Light curve of the 2003 outburst of XTE J1814—338. 
Legends are same as in Fig. ^ 


HETEJ1900.1-2455 



Time (MJD) 


SAX J1808.4-3658 (1998 Outburst) 




Time (MJD) 


Figure 8. Light curve of the 1998 and 2002 outbursts of 
SAX J1808.4-3658. Legends are same as in Fig.[^ 


Figure 7. Light curve of HETE J1900.1-2455 from June 2005 to 
January 2007. Legends are same as in Fig.[^ Note that there are 
additional pulse detections beyond ~ 54000, however, these are 
tentative | |Patrun^|2012[ | and do not change the measured flux 
range. 


The distance to the source is estimated at < 8 kpc, by 
assuming the measured luminosity during the type I X-ray 
burst is Eddington limited ( |Strohmayer et aL]|2003| ). The 
resulting magnetic field range is 1.6 x 10' G < B < 7.8 x 
10® G. 


3.7 HETE J1900.1-2455 


HETE J1900.1-2455 was discovered through a bright type 
I X-ray burst observed with the High Energy Transient Ex¬ 


plorer 2 (HETE-2) in June 2005 (Vanderspek et al. 20051 
and 377 Hz pulsations were observed with RXTE quickly 
thereafter ( Morgan et al.p005 Kaaret et al.|[2006 |. Unlike 
the other AMXPs, which show outbursts that last for weeks 
to months, HETE J1900.1-2455 has been active since discov¬ 
ery and is yet to return to quiescence. Persistent pulsations 
have been reported to occur during the first 20 days of the 
outburst, after which only intermittent pulsations have been 
seen ( Patruno|2012 1 . 

We have analysed all archival RXTE observations 


of this source and find that the observations of highest 
and lowest flux with pulsations are on MJD 53559.5 and 
MJD 53573.8, respectively. 

Since the source has shown continuous activity since 
its discovery, there is no observation during quiescence from 
which the background flux can be estimated. However, as 
the source is well away from the galactic centre (galactic 
coordinates: I = 11.3°, b = —12.9°), the the background firrx 
is expected to be comparatively small. We therefore neglect 
the background contribution for this source, noting again 
that this leads us to slightly overestimate the source flux, 
which yields a less stringent constraint on pmax- 

The distance to the source is taken to be d ~ 5 kpc, 
based on photospheric radius expansion of a type I X-ray 
burst ( Kawai fc Suzuki|2005 1. The magnetic field range ob¬ 
tained isTT63n'o'^~G~<'B'< 1.0 X 10® G. 


3.8 SAX J1808.4-3658 


SAX J1808.4-3658 was discovered with BeppoSax in in 1996 
( in ’t Zand et al)]|1998 ) and the detection of 401 Hz pulsa¬ 
tions in 1998 with RXTE made it the first known AMXP 
| |Wijnands fc van der Klis||199^ . SAX J1808.4-3658 has 
been observed in outbursts with RXTE six times. 

We find the highest flux with pulsations to be in the 
2002 outburst on MJD 52563.2 and the lowest flux with 
pulsations in the 1998 outburst on MJD 50936.8. The low- 
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IGRJ17498-2921 



Time (MJD) 


Figure 9. Light curve of the 2011 outburst of IGR J17498-2921. 
Legends are same as in Fig. ^ 


XTEJ1751-305 



Time (MJD) 

Figure 10. Light curve of the 2002 outburst of XTE J751-305. 
Legends are same as in Fig. ^ 


est flux with pulsations is observed toward the end of the 
1998 outburst (Fig.j^. Although the flux of this observation 
is only slightly higher than that of the background observa¬ 
tions, the pulse detection is very significant and yields a 
phase that is consistent with the timing model of this out¬ 
burst. If pulsations are still present at the same amplitude 
in the other low flux observations, we would not be able to 
detect them due to the low count rate. We therefore consider 
the low flux observation with pulsations to be background 
dominated. 

The distance to SAX J1808.4-3658 is d = 3.5 kpc, and 
was derived from photospheric radius expansion in type I X- 


ray bursts Galloway & Gumming (20061. With this distance 
we obtain a magnetic field range of 1.4 x 10^ G < B < 

1.8 X 10® G. 


3.9 IGR J1T498-2921 

IGR J17498-2921 was discovered in August 2011 with IN¬ 
TEGRAL dGibaud et aL][20TT|, following which pulsations 
at 401 Hz were reported by Papitto et al. (2011bI. The light 
curve of this outburst is shown in Fig. |9| 

Since the source is close to the galactic centre (galac¬ 
tic coordinates: I = 0.16°, fe = ~1°), the observations have 
a large X-ray background contamination, with the lowest 
observed flux with pulsations again background dominated. 

The distance to IGR J17498-2921 is estimated at ~ 
8 kpc, based on photospheric radius expansion during a type 
I X-ray burst ( Falanga et al.|2012 K We then find a magnetic 
field range of TO x 10' G < B < 1.6 x 10® G. 


We estimate the background from the last RXTE obser¬ 
vation on MJD 52380.7, when the pulsations were no longer 
detected (see Fig. [To|). 

The distance to this source is not well-defined. IMark-l 


wardt et al. 


(2002a I constrain the distance to U 7 kpc by 
equating predicted mass transfer rates | |Rappaport et ah] 
1983 King et al. 19971 to values inferred from X-ray ob¬ 


servations. Papitto et al. (20081 instead compare the spin- 


frequency derivative to models for angular momentum ex¬ 


change (Rappaport et al. 20041, and so constrain the dis¬ 


tance to 6.7 — 8.5 kpc. For our analysis, we take a dis¬ 
tance of 7 kpc, which results in a magnetic field range of 
2.5 X 10^ G < B < 1.1 X 10® G. 


3.11 SAX J1748.9-2021 

SAX J1748.9-2021 is located in the globular cluster 
NGC 6440 (in ’t Zand et al. 19991, and was observed in 


outburst by RXTE in 1998, 2001, 2005 and 2010 ( 

in’t Zand 

et al.|19991|in’t Zand et al.|2001t|Markwardt & Swank|20051 

Patruno et al. 

2010b 1 . Among these outbursts 442 Hz in- 


termittent pulsations have been detected in 2001, 2005 and 


2010 (Gavriil et al. 2007 Altamirano et al. 2008 Patruno 
et al.|2010b K 


We find that the highest flux with pulsations occurs dur¬ 
ing the 2010 outburst on MJD 55222.5, while the lowest flux 
with pulsations is seen in the 2001 outburst on MJD 52198.3 
(light curves shown in Fig. EH- 

As the source is associated with a globular cluster, 
its distance is comparatively well constrained to 8.2 kpc 
(| Valenti et al.||2007|). We find a magnetic field range of 
T9 X 10" G < B < 3.8 x 10® G. 


3.10 XTE J1751-305 

XTE J1751-305 and its 435 Hz pulsations were discov¬ 


ered with RXTE in April 2002 (Markwardt et al. 2002aI. 


A brief second outburst was detected in 2009 ( Markwardt 
|et al.|[^09p , which was coincident with on-going activity of 
IGR J17511-3057 in the same field of view. Since RXTE is 
not an imaging detector, the flux contribution of these two 
sources cannot be separated, so we restricted our analysis to 
the 2002 outburst of XTE J1751-305 only. 


3.12 Swift J1749.4-2807 


Swift J1749.4-2807 was first detected in June 2006 (Schady 


et al. 20061, but its 518 Hz pulsations were not found un¬ 


til the second outburst in April 2010 (Altamirano et al. 


2010c 2011 We find both the highest and lowest flux with 


pulsations occur in this second outburst on MJD 55300.9 
and 55306.7, respectively. 

The source is located close to the galactic centre and 
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SAX J1748.9-2021 (2001 Outburst) 



55220 55230 55240 55250 


Time (MJD) 


Figure 11. Light curves of the 2001 and 2010 outbursts of SAX 
J1748.9-2021. Legends are same as in Fig.[2 
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Figure 12. Light curve of the 2010 outburst of Swift J1749.4- 
2807. Legends are same as in Fig. 


has a strong contaminating X-ray background flux. The 
absence of pulsations (Fig. 121 indicate that the source 
was no longer accreting for the last three RXTE observa¬ 
tions (> MJD 55307), which is confirmed by a source non¬ 


detection with both Swift and INTEGRAL (Ferrigno et al. 
|2011[ ). We therefore use the last RXTE observation to mea¬ 
sure the background flux. As this background flux is compa- 


IGR J00291+5934 



Figure 13. Upper panel: Light curve of the 2004 outburst of 
IGR J00291+5934. Legends are same as in Fig. ^ 


rable to the lowest observed flux with pulsations, we cannot 
confidently estimate the source contribution and to be con¬ 
servative we again calculate the upper limit on the magnetic 
held without using the background. 

The distance to the source is d = 6.7 ± 1.3 kpc, which 
was inferred from the luminosity of a suspected type I X-ray 
burst i Wijnands et al.|2009 l. Adopting the central value for 
the distance and the measured hux we obtain a magnetic 
held range of 1.1 x 10^ G < B < 7.7 x 10® G. 


3.13 Aql X-1 


In 18 outbursts across ~ 15 years, Aql X-1 has shown its 
550 Hz pulsations only for a single episode of about 150 sec¬ 


onds (Casella et al. 2008 Messenger & Patruno 2014 1 . Hence 


we cannot measure a hux range for the presence of pulsations 
for this source. Instead we calculate both lower and upper 
limits on the magnetic held from the same measured hux, 
thus obtaining a very conservative estimate of the allowed 
magnetic held range. 

The distance to the source is 4.4 — 5.9 kpc and was 
obtained from a photospheric radius expansion during a type 
I X-ray ( Jonker fc Nelemans||2004 1. For a distance of 5 kpc 
and the measured hux we obtain a magnetic held range of 
4.4 X 10'^ G < H < 3.1 X 10® G. 


3.14 IGR J00291-I-5934 

IGR J00291-I-5934 was discovered with INTEGRAL in De¬ 
cember 2004 (Eckert et al. 2004 Shaw et al.||2005 l and its 
599 Hz pulsations were detected in follow-up RXTE obser¬ 
vations ( [Markwardt et ar]|2004[ ). Outbursts were detected 
again in August and September 2008 (Chakrabarty et al. 
2008 Lewis et al.||2008 l. We hnd both the highest and the 
lowest hux with pulsations occur in the 2004 outburst. 


It is evident from the light curve (Fig. 131 that the 


source gradually decays to the background level, which is 
due to an intermediate polar in the same held of view 
( Falanga et al.p005b ). Because the lowest hux with pulsa¬ 
tions is again comparable to the estimated background hux, 
we assume the low hux observation is background dominated 
to get a more conservative magnetic held estimate. 
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1 % Eddington rate - - - 
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Figure 14. Magnetic field vs spin frequency for AMXPs anal¬ 
ysed in this work. The grey lines show the range of estimated field 
strengths of the sources from Table. T he black lines represent 
the revised upper limits following eq. |ll| which gives tighter con¬ 
straints of the field strength estimates. The dashed and dashed- 
dotted curves show the equilibrium spin rate for two different 
accretion rates. 


4.2 Masses and radii of AMXPs 


Currently there are no reliable estimates available for the 
masses or radii of AMXPs. In calculating the field strengths 
presented in Table we assumed the neutron stars to be of 
canonical mass M ~ IAMq and radius ~ 10 km. However, 


theoretical calculations of mass-radii relations (e.g. Lattimer 
fc Prakash||2001 Lattimei]|2014 1 predict that both parame¬ 
ters may vary over a wide range of M ~ 0.2 — 2.9Mq and 
7? ~ 10 — 15 km. 

For a more massive neutron star of ~ 2Mq the up¬ 
per limit on the magnetic dipole moment will increase by 
~ 12.6% (/imax oc whereas the lower limit will de¬ 
crease by 8.2% (/imin oc from the values reported in 

Table resulting in a broadening of the estimated range. 

For a neutron star of larger radius, e.g. ~ 15 km, the 
upper limit on the magnetic dipole moment will scale as 
Mmax oc B}^^, causing an increase of 22.5% with respect 
to the values we report. The lower limit to the magnetic 
dipole moment scales with radius as /rmtn oc such that 

assuming a 15 km radius increases the lower limit by 149%. 
For a 7? = 15 km neutron star we then find that the allowed 
range of magnetic field strength moves to higher values and 
covers a slightly broader rang^ 


The distance to IGR J00291-I-5934 has been esti¬ 
mated in several ways. From the long term average accre¬ 
tion rate, Galloway et al. (2005) constrain the distance to 
^ 4 kpc. Gomparing the observed quiescent flux to that of 
SAX J1808.4-3658, [Jonker et al.| ( |2005| ) estimate the dis¬ 
tance to be 2 — 3.6 kpc. Similar estimates were also obtained 
by Torres et al. (20081, who report d — 1.8 — 3.8 kpc by 
modeling the light curve. In this work we adopt the central 
distance of d = 3 kpc, which gives a magnetic field range of 
8.5 X 10® G < B < 1.9 X 10® G. 


4 DISCUSSION 

In this section we discuss the various sources of uncertainty 
involved in the magnetic field estimates presented in this 
work. 


4.1 Distance estimates 

The largest uncertainties in our field estimates are due to 
poor constraints on the distance. Errors in the distance af¬ 
fect the luminosity estimates and cause a systematic uncer¬ 
tainty that scales both the upper and lower limit by the 
same factor (T^min/max CK d). The sources associated with 
globular clusters are not greatly affected by this uncertainty 
as they have well measured distances. Some AMXPs exhibit 
thermonuclear bursts with photospheric radius expansion, 
which allows for tight constraints on the distance. For the 
remaining sources, however the distance measurements are 
less certain and the resulting systematic uncertainty has a 
more prominent role. In Table|^we explicitly report the dis¬ 
tances used in calculating the upper and lower limits, such 
that our estimates can be easily adjusted should a more ac¬ 
curate measure distance be obtained for one of these sources. 


4.3 Disk-magnetosphere interaction 


In this paper we made the assumption that the neutron star 
magnetic field is dipolar in nature, which may not be valid. 
Near the truncation radius the pressure the disk exerts on 
the magnetosphere can affect the field geometry, and near 
the neutron star surface this may be further complicated 
by higher multipole moments that could dominate over the 
dipole component. We parametrized the effect of the field 
shape with the factor ys (see eq. |^, which we assumed 
to vary over a wide range of 0.01 — 1 to account for such 
uncertainties. This range is consistent with the results of 
numerical simulations, which suggest the largest range of 
truncation radii is rt ~ (0.5—1.2)rA (Zanni &; Ferreira|2009 


Romanova et al.]|2008l, implying ys — 0.06 — 1.3. 


• Recent MHD simulations (Zanni & Ferreira 2013 
Kulkarni fc Romanova|[2013| [Long et al.||2008[ |R ornanov^ 


et al. 


20081 have confirmed that complex field topologies 


arise at the inner edge of the disk. [Kulkarni fc Romanova| 
(20131 show that if the truncation radius is in the range 


2.577s < rt < 577s, the non-dipolar field structure results in 
a modified expression of the truncation radius 


n = 1.0677, 


\gmm^ri) 


( 




1/10 


( 11 ) 


This modification is relevant when estimating the upper 
limit to the magnetic field, for which we take the disk to 
be truncated at the co-rotation radius, which normally falls 
within the range of radii where eq. is applicable. Using 


^ The field strength range covers more than an order of magni¬ 
tude, so despite the large fractional change, the absolute shift of 
the lower limit is smaller than that of the upper limit. 
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this expression for the truncation radius we obtain a modi- 
hed boundary layer parameter as 


’•'"■““(lorH;) (l 


M 


4 Mq 


- 1/2 


Ra 


10 km 


\ 3/2 

j (12) 


which, instead of conservatively assuming 0.01, may be used 
in eq. to obtain a more constraining upper limit on the 
magnetic dipole moment. 

• Near the neutron star surface higher multipole moments 
of the magnetic held may be stronger than the dipole com¬ 


ponent. As shown in numerical studies (Long et al. 2007 


20081 such complex non-dipolar held conhgurations strongly 


affect the inner accretion geometry, but the ehect on ys is 
not well established. If further theoretical considerations can 
constrain ys to a smaller value than our conservative as¬ 
sumption of ys = 1, then this could be used to tighten the 
lower limit on the magnetic held strength {Bn 


-l/2x 

« 7s )• 


In Fig. [^we plot the conservative ranges of the sur¬ 
face held strength following the analysis in Sec. |2.1| in grey 
and present the more constraining estimates based on eq. 
in black. The dashed lines represent the equilibrium spin 
(Alfven radius at co-rotation) for a mean long term accre¬ 
tion rate 

-6/'^ / R \ -18/'^ 


Van = 441 Hz ( —^ 
eq ^10 G 


x(-) 

KmeJ 


3/7 


M 


1.4 Mq 


Rs 

10 km 

5/7 


(13) 


where Ms = 1.5xl0 * Mq yr ^ is the Eddington accretion 
rate. 

• In our analysis we assumed that channelled accretion 
onto the neutron star can only take place when the disk is 
truncated inside the co-rotation radius. However, the large 
range of X-ray luminosities, and accordingly mass accretion 
rates, observed for AMXPs suggests that mass accretion 
onto the neutron star might persist even when the inner 
edge of the disk moves outside the co-rotation radius (see 
e.g. Rappaport et al.|2004 l, which indeed appears to be con- 
hrmed by observation | Bult & van der Klis 


As pointed out by Spruit & Taam (19931 the inner edge 


20151. 


of the disk must have receded to rt ~ 1.3 before the cen¬ 
trifugal force is strong enough to accelerate matter beyond 
the escape velocity and thus drive an outflow. If we consider 
the possibility that accretion may still occur for radii up-to 
1.3 Cc, then we find that the upper limit to the magnetic held 
strength increases by 58% (Bmax oc {rtjrcfi at most. The 
lower limit, being independent of the co-rotation radius, is 
unaffected. 


4.4 Observational sampling 

• In order to determine the upper limit to the magnetic 
held strength we consider the lowest observed hux for which 
pulsations are signihcantly detected. However, as the hux de¬ 
cays, the signal-to-noise ratio also decreases, such that the 
non-detection of pulsations could be due to limited statis¬ 
tics. That is, the pulsations may persist below our detection 
limit. This concern is particularly relevant to XTE J1807.4- 
294, IGR J17498-2921, and IGR J17511-3057 in which pul¬ 
sations are detected at approximately the same hux as the 


background level. If pulsations are still present at a very 
low level, and accretion is ongoing at even lower luminosi¬ 
ties, then our estimates for these sources are overly con¬ 
servative as we are overestimating the upper limit to the 
magnetic held strength. Since Hmax oc ^ future de¬ 

tection of pulsations at a lower luminosity than reported 
in this work can therefore be used to further constrain the 
range of magnetic held strengths. This is especially relevant 
for other X-ray satellites that have better sensitivity and 
a lower background contamination such as XMM-Newton, 
or for a future ASTROSA T ( |Agrawal|[200^ |Singh|[20l4| or 
LOFT ( Feroci et al. [2012 1 mission. 

• For the lower limit to the magnetic held strength we 
use the highest observed hux for which pulsations are de¬ 
tected. As we note in Sec. the pulsations are expected to 
disappear at high hux when the disk extends to the neu¬ 
tron star surface, yet this is never observed. For all sources 
the highest observed hux considered always shows pulsa¬ 
tions, which implies that the inner edge of the disk never 
extends down to the neutron star surface. Indeed some of 
the better sampled AMXPs (e.g. SAX J1808.4-3658 and 
IGR J00291-I-5934) show peak luminosities that vary by a 
factor of 2 between outbursts. If a future outburst reaches a 
higher peak luminosity than considered in this work, it will 
increase the lower limit as Bmin oc and thus further 

constrain the allowed range of magnetic held strength. 


4.5 Luminosity estimates 

• To calculate the hux we consider the 3-20 keV X-ray 
band, which we convert to the bolometric hux by applying 
a correction factor (Gilfanov et al.|1998 Casella et al.|2008 


Galloway et al.|2002 Campana et al.|2003 1. To be conserva¬ 

tive in our estimates of the held strength we used a Cboi = 2 
for the upper limit and Cboi = 1 for the lower limit. How¬ 
ever, if the correction factor is well constrained then this 
approach is overly pessimistic. For many accretion powered 
pulsars the correction factor tends to be within ~ 10% of 2 
( Galloway et al.|2008 l, such that the error introduced in the 
magnetic held estimate is only ~ 2% (Bmax oc (eboi/2)^^^, 
Bmin oc (cboi/1)^^^). If we adopt the same bolometric cor¬ 
rection factor for the lower limit also, it we hnd the our 
estimates can improve by up-to < 30%. 

• Another source of uncertainty comes from the back¬ 
ground contribution to the measured hux. We estimate the 
background contribution from the RATE observations at the 
end of an outburst, when the source no longer shows pulsa¬ 
tions and has presumably returned to quiescence. If this es¬ 
timate contains residual source emission, or the background 
contribution itself is variable, then this approach introduces 
an error in our held strength limits. 

Because the highest observed hux with pulsations is al¬ 
ways much higher than the background the lower limit to 
the magnetic held strength will not be greatly affected by 
the uncertainty in the background estimate. The lowest ob¬ 
served hux with pulsations, however, is often comparable to 
the background contribution, so the effect on the upper limit 
needs to be considered carefully. 

For some sources the background hux is sufficiently lower 
than the minimum pulsating hux (as shown in Table 
that the effect of the background correction amounts to 
only a small change in the held estimates (e.g. ~ 8% 
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for XTE J1814-38, ~ 12% for Swift J1756.9-2508). For 
sources where background flux could not be measured 
(HETE J1900.1-2455 and XTE J0929-314) the error intro¬ 
duced by neglecting the background is expected to similarly 
be only a few percent. 

The remaining sources have a comparatively high back¬ 
ground contribution, such that the estimated background 
flux is a large fraction of the minimum pulsating flux. For 
these AMXPs we conservatively took the lowest observed 
flux with pulsations as an upper limit to the real flux at 
which pulsations and thus accretion stops. We then calcu¬ 
lated the upper limit on the magnetic field strength with¬ 
out adjusting for the background. Accounting for the back¬ 
ground would lower the source flux estimate by up-to a 
factor of two, and thus improve our estimates by roughly 
< 40%. Further improvement might be achievable with more 
sensitive instrumentation (as noted in the previous section). 

• To convert the observed flux to the source luminosity 
we assumed an isotropic emission process, however, the flux 
includes contributions from the hotspot, which may have a 


significant anisotropy (Poutanen & Beloborodov 2006 Vi- 


[ironen &: Poutanen||2004[ [Poutanen fc Gierlihski]2003[ ). 

The effect of anisotropic emission is not at all clear. The 
degree of anisotropy of the hotspot emission depends on as¬ 
sumptions of the emission process and can vary by a factor 
of ~ 2 ( Poutanen fc Beloborodov|2006 1. Furthermore, what 
fraction of the total flux is affected by this will depend on 
the size, shape and position of the hotspot and is subject 
to considerable uncertainty. At best this effect applies only 
to the pulsed component of the emission (~ 10%) and thus 
introduces a systematic error in our estimates of only a few 
percent. At worst most of the observed flux originates from 
a large hotspot that has a slightly beamed emission pattern. 
In that case the allowed magnetic field range may show a 
systematic shift of up-to ~ 40%. 


5 COMPARISON WITH PREVIOUS WORKS 

In this work we estimated upper and lower limits to the 
magnetic field strength of all accreting millisecond X-ray 
pulsars (AMXPs) observed with RXTE. We assume that the 
detection of X-ray pulsations signifies ongoing magnetically 
channelled accretion. Thus, by associating the range of lumi¬ 
nosity for which pulsations are detected with the expected 
extent of the disk truncation radius, we have constrained the 
magnetic dipole moment of the neutron star. The obtained 
equatorial surface magnetic field strengths of the 14 AMXPs 
analysed are presented in Table 

Our magnetic field strength estimates are subject to a 
number of uncertainties, which were discussed in the previ¬ 
ous section. We note that we have chosen most of the uncer¬ 
tain parameters such that we obtained a conservative range 
for the magnetic moment. If parameters such as ys or Cboi 
can be established more accurately, they will tighten the 
constraints further. Errors in other parameters, such as the 
distance to the source, introduce a systematic shift in our 
results. Refinement in the measurement of these parameters 
would affect both the upper and the lower limits the same 
way. Similarly, the uncertainties in the estimation of the 
background fluxes also affect both limits to the magnetic 
field. If better estimates are available from more sensitive 


instruments and also if there are future outbursts with flux 
ranges wider than those in the current work, the constraints 
on the magnetic field in Table can be easily updated by 
correcting the fluxes presented in Table and recomputing 
the limits to the dipole moment using eq7[^ and eq. |10| 

For 5 of the 14 considered AMXPs (Swift J1749.4-2807, 
IGR J17511-3057, NGC 6440 X-2, XTE J1807.4-294 and 
IGR J17498-2921) we obtain constraints on the magnetic 
field strength for the first time. For the other 9 AMXPs field 
strength estimates have been previously reported. Below we 
discuss some of the techniques used to obtain those estimates 
and how they compare with limits we report here. 


(i) Vacuum dipole radiation in quiescence: For some 
AMXPs the spin frequency has been measured for successive 
outbursts that are months or years apart. Measurements of 
spin-down during the intervening periods of quiescence, can 
then be used to estimate the magnetic field strength by as¬ 
suming the spin-down is due to magnetic dipole emission. 
Such estimates have been obtained for IGR J00291-I-5934 
||Patruno|2010[ [Hartman et al.|201ip , XTE J1751-305 (jRig-j 
gio et al.pOlla l, SAX J1808.4-3658 ( Hartman et al.||2008 


2009 


Patruno et al.|20l2 l and Swift J1756.9-2508 (Patruno 


20101. While the radio emission associated with magnetic 


braking has not yet been detected in any of these sources, the 
recent discovery of millisecond neutron stars pulsating alter¬ 
nately in X-rays and radio (jPapitto et al.|2013b| [Archibald] 


et al.|2014 Papitto et al.|20T5 1 provides some evidence that 

this interpretation of measured spin-down is correct. 

Our magnetic field strength estimates appear to be sys¬ 
tematically lower than those obtained through quiescent 
spin down, although we note that given the systematic un¬ 
certainties discussed in the previous section, the results of 
these two approaches are roughly consistent. 

(ii) Quiescent luminosity estimates: For some LMXBs e.g. 


SAX J1808.4-3658, Aql X-1 (Di Salvo & Burderi 2003b|, 
KS 1731-260 ([Burderi et al.|2002| and XTE J0929-314 ([Wl 


jnands et al. 20051, limits to the magnetic field strength have 


been inferred from measurements of the quiescent X-ray lu¬ 
minosity. However, given the very low count rates in the 
quiescent phase, and a poor understanding of which physical 
mechanism governs the radiation process, these methods of¬ 
fer less reliable constraints on the dipole moment compared 
to other approaches. While for Aql X-1 no other independent 
confirmation of the upper limits exists, the upper limits for 
SAX J1808.4-3658 and XTE J0929-314 obtained through 
quiescent luminosity methods are comparable to the upper 
limits we report here. 

(iii) Accretion induced spin down estimates: For some sys¬ 
tems the magnetic field strength has been estimated by com¬ 
paring observed rate of spin down during ongoing accretion 
to theoretical estimates of magnetic torque. For example, for 


XTE J1814-338 Papitto et al. (20071 assume that observed 


pulse frequency variations are caused by spin down due to 
the torque applied by an accretion disk that is truncated 
near the co-rotation radius. Following theoretical calcula¬ 


tions of accretion induced spin down torques ( Rappaport 
[et al.|2004[ ), the authors estimate the surface magnetic field 
to be ~ 8 X 10® G, which is comparable our conservative 
upper limit to the field, and significantly higher than the 
upper limit we obtain using eq. [12[ However, given their 
simplifying assumptions regarding the magnetic field topol- 
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ogy at the Alfven radius, and the considerable uncertainty 
in interpreting pulse frequency variations as spin variations 
(Patruno et al. 2009el, accretion induced spin down esti¬ 
mates are less robust than those we obtain. 

(iv) Burst oscillations: In some accreting pulsars the 
phase of burst oscillations in type I X-ray bursts is locked to 
the phase of the accretion powered pulsations (XTE J1814- 
338 [Watts et al.]|2008| IGR J17480-2446 [Cavecchi et al.| 


2011 ). 


Cavecchi et al. 


(20111 argue that if this phase-locking 


is due to magnetic confinement of the flame propagation 
front, then it would require a field strength of 5 x 10® G. 
Even given the systematic uncertainties that enter in our 
estimates, such a large magnetic field strength would be dif- 
hcult to reconcile with the lower upper limit we obtain for 
XTE J1814-338. 

(v) Spectral state transitions in LMXBs: In another ap¬ 
proach, spectral state transitions have been used to iden¬ 


tify the onset of the propeller regime (Matsuoka & Asai 
2013 Asai et al.|2013 1. These authors argue that a spectral 


change and a fast decline of luminosity towards the end of 
an outburst indicates the onset of the propeller regime and 
thus an accretion disk that is truncated at the co-rotation 
radius. Such spectral state transitions have been investi¬ 
gated for only a handful of LMXBs (e.g. Aql X-1, 4U 1608- 
52, XTE J1701-462), but give limits on the magnetic held 
strength of Aql X-1 that are tighter than those we obtained. 
However, there is no clear evidence that these transitions are 
indeed caused by a propeller effect. In fact, the observation 
of similar state transitions in black hole binaries seems to 
suggest otherwise (Jonker et al.|2004). 


To conclude, the magnetic held estimates we obtain 
agree with most of the other indirect methods to within 
an order of magnitude. The large uncertainties on many pa¬ 
rameters, as well as the uncertainty in the underlying as¬ 
sumptions, introduce a signihcant spread in the range of 
held strengths inferred via each approach, with no single 
technique being more robust than the others. Nonetheless, 
all alternative methods discussed in this section require that 
the AMXP be observed during a specihc state of its outburst 
(e.g during quiescence or a spectral state transition). Since 
observations of such special states are not available, or may 
not even exist, for all sources, these methods are not suit¬ 
able for studying the population. By contrast, the method 
used in this work applies to any AMXP, making it more reli¬ 
able for comparing the held strengths of the population and 
understanding the evolutionary processes which lead to the 
formation of AMXPs. 
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Table Al. Adopted Tasc values for the outbursts of NGC 6440 
X-2. 


NGC 6440 X-2 


Outburst 

Tasc 

Outburst 1 

55042.817 

Outburst 2 

55073.034 

Outburst 3 

55106.012 

Outburst 4 

55132.907 

Outburst 5 

55276.625 

Outburst 6 

55359.470 

Outburst 7 

55473.854 

Outburst 8 

55584.714 

Outburst 9 

55871.231 


A6 XTE J1814-338: 


The RXTE data of XTE J1814-338 is given by program IDs 
P80(138/145/418) and P92054. The optical position used to 


barycenter the data is given by Krauss et al. 


timing solution is taken from Papitto et al. 


(2005a I and the 


(20071. 


A7 HETE J1900.1-2455: 

The source HETE J1900.1-2455 has a long history of ac¬ 
tivity and thus its RXTE data is spread over many pro¬ 
gram IDs: P910(15/57/59), P91432, P92049, P93(030/451), 
P940(28/30), P95030 and P96030. For source coordinates we 


use the optical position of Fox (20051. The timing solution 
is given by [Patruno (20121. 


APPENDIX A: DETAILS OF THE TIMING 
ANALYSIS 

Al Swift J1T56.9-2508: 

The two outbursts of Swift J1756.9-2508 were observed 
with RXTE under program IDs P92050 & P93065 (2007) 


and P94065 (2009). We barycenter this data using the best 

known coordinates of Krimm et al. 

(2007bI. For the timing 

solution we adopt the ephemeris of 

Patruno et al. (2010a 1. 


A2 XTE J0929-314: 

We correct the RXTE data (program ID P70096) to the 
Solar system barycenter using the Chandra source position 
of Juett et al. (20031. The data was folded using the timing 
solution of lacolina et al. (20091. 


A3 XTE J1807.4-294: 


The RXTE data of the outburst of XTE J1807.4-294 is 
given by program IDs P70134 and P80(145/419). We use the 
Chandra coordinates of Markwardt et al. ( |2003 1 to barycen¬ 
ter the data. We use the system ephemeris of Riggio et al.| 


(20081 to fold the data. 


A4 NGC 6440 X-2: 

The nine short outbursts of NGC 6440 X-2 that were 
observed with RXTE are given by program IDs P94044, 
P94315 and P95040 ( Patruno fc D’Angelo|2013 1. We correct 
the data to the Solar System barycenter using the Chandra 


position of Heinke et al. (20101. The timing solution for the 


hrst outburst is given by Altamirano et al. (2010al. This 


timing solution is also used the later outbursts, but with 
the locally optimized time of ascending node values given in 
Table lAll 


A5 IGR J17511-3057: 

The RATE data of the 2009 outburst of IGR J17511-3057 is 
given by program ID P940(41/42). We barycenter the data 
using the Chandra position given by Nowak et al. ( 2009[|. 
The timing solution of this source is given by | Riggio et al.| 
(2011b|). 


A8 SAX J1808.4-3658: 


The RXTE program IDs for SAX J1808.4-3658 are given 
by P30411 (1998), P40035 (2000), P70080 (2002), P91056 & 
P91418 (2005), P93027 & P93417 (2008) and P96027 (2011). 
When considering the outbursts of SAX J1808.4-3658 we 
exclude the prolonged outburst tail which represents an un¬ 
usual disk state ( [Patruno et al.|[2009d[ |2015| ), so the out¬ 
burst of 2000 is entirely omitted from the analysis. Source 
coordinates for barycentering are taken from [Hartman et al. 


([2008p. For the ephemeris we adopt the solutions of 


et al. (2008 20091 and Patruno et al. (2012| for the respec- 


Hartman 


tive outbursts of 1998, 2002 and 2005; 2008; and 2011. 


A9 IGR J17498-2921: 

RXTE data for IGR J17498-2921 is given by program ID 
P96435. For the source position we adopt the Chandra co¬ 
ordinates of Chakrabarty et al. (20111. The timing solution 
is given by Papitto et al. (2011a|. 


AlO XTE J1751-305: 


We only analyze the RXTE data of the 2002 outburst (see 
main text for details) which is given by programs P70131 
and P70134. To barycenter the data we use the Chandra 
position of Markwardt & Dobrzycki ( 2002[). For the orbital 
ephemeris we use the timing solution of [Markwardt et al.| 
(2002b|). 


All SAX J1748.9-2021: 

The RXTE data of SAX J1748.9-2021 is given by pro¬ 
grams P30425, P60035 & P60084 (2001), P91050 (2005) and 


P94315 (2010). We use the Chandra position of [in’t Zand 
et al. (200ip for barycentering the data. The timing solution 


of Patruno et al.[ (2009a i was used to fold the data of the 


1998, 2001 and 2005 outbursts and that of [Patruno et ah] 
(2010b I for the 2010 outburst. 


A12 Swift J1749.4-2807: 


We use the data of RXTE program P95085. For barycen¬ 
tering we adopt the X-ray coordinates of [Wijnands et ah] 
( [2009p, obtained wit h XMM-Newton. We use the ephemeris 
of Altamirano et al. (20111 to fold the data. 
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A13 IGR J00291+5934: 

We use RXTE data from programs P90052 and P90425 for 
the 2004 outburst and program P93013 and P93435 for the 
2008 outburst. We use radio coordinates of |Rupen et al. 


(20041 to barycenter the data. The timing solution of 


truno 


Pa- 


(20101 was used to fold the data. 


HF spectrum was fit with a power-law continuum and ther¬ 
mal blackbody component whereas the LF and BE spectra 
were fit with a simple power-law model. The hydrogen col- 


umn density was fixed to Nh = 5.6 x 10^^ cm ^ (Falanga 


et al.||2005al. The resulting fit parameters are presented in 
TablelB3l 


APPENDIX B: DETAILS OF THE SPECTRAL 
ANALYSIS 

For the spectral analysis we extracted the data of all Xenon 
layers of PCU2, and for observations with a poor signal to 
noise ratio we combined data from all active PCUs. Due 
to poor signal to noise some spectral fits at low flux have 
a reduced-y^ of much less than 1. As our work does not 
focus on obtaining the most accurate spectral model, but 
rather on measuring the flux, our final results will not be 
signihcantly affected by the model uncertainty of these fits. 

Our analysis of the RXTE/PCA data was performed for 
the 3 — 20 keV energy range using XSPEC version 12.7.1. 
The galactic absorption was modelled with the TBABS 


model (Wilms et al. 20001. Since RXTE instruments can¬ 


not properly constrain the galactic absorption due to neu¬ 
tral hydrogen in the lower energy range, we fix Nh to values 
obtained from literature. The errors quoted are for a 95% 
confidence limit. 


B1 Swift J1756.9-2508: 

For ObsIDs 94065-02-01-05 (high flux with pulsations, HF) 
and 94065-06-02-03 (low flux with pulsations, LF) data were 
extracted from PCU2 whereas for 94065-06-03-02 (back¬ 
ground estimate, BE), data from PCU 0, 2 and 4 were com¬ 
bined for the spectral analysis. The HF and BE spectra were 


fit with a simple absorbed power-law model XSPEC (Ar- 
naud|1996 1. The LF spectrum was fit with a power-law con¬ 
tinuum and a thermal black body component. The hydrogen 
column density accounting for galactic absorption was fixed 
to Nh = 5.4 X 10^° cm“^ (Krimm et al. 2007b I. A Gaus¬ 


sian emission component centred at ~ 6.5 keV (width fixed 
to 10“^ keV) was required to obtain statistically acceptable 
hts. The results of the spectral analysis are presented in Ta¬ 
ble [SI 


B4 NGC 6440 X-2: 

For ObsIDs 94044-04-02-00 and 96326-01-36-00 data was ex¬ 
tracted from PCU2 whereas for ObsID 96326-01-40-01 data 
was taken from PCU 1 and 2. The h ydrogen colu mn density 
was fixed to Nh = 5.9 x 10^^ cm~^ ( Harris|l996 K The data 
were fit with a simple absorbed power-law for all ObsIDs. 
The fit parameters are presented in Table |B^ 


B5 IGR J17511-3057: 

For ObsID 94041-01-01-02 (HF) the data were extracted 
from PCU 2. For ObsIDs 94042-01-03-04 (LF) and 94042- 
01-02-05 (BE) spectral analysis was performed by combining 
data from PCUs 2, 3 and PCUs 2, 4 respectively. The hy¬ 
drogen column density was fixed to Nh = 1 x 10^^ cm“^ 


(Papitto et al. 2010[ Paizis et al. 20121. The HF was well 
described with a power-law continuum and thermal black- 
body component, whereas the LF and BE spectra were fit 
with a simple absorbed power law. A Gaussian feature at 
~ 6.5 keV was added to the HF and BE spectra to improve 
the fits. The resulting parameters are presented in Table |B^ 


B6 XTE J1814-338: 

Since the observed count rate of XTE J1814-338 is small 
(< 40ct s“^ for the entire outburst), data from all available 
PCUs were combined for all spectra. For ObsID 80418-01- 
03-08 (HF) the data were extracted from PCU 0, 2 and 
3. For 80418-01-07-08 (LF) data from PCUs 0, 1, 2 and 3 
and for 80418-01-09-00 (BE) PCUs 0, 1 and 2 were used. 
The hydrogen column density was fixed to Nh = 1.67 x 


lO^'^ cm ^ ( Krauss et al.|[^05b l. 


The HF spectrum was fit with a power-law continuum 
and a blackbody component, whereas the LF and BE spectra 
were fit with a simple absorbed power-law model. The details 
of the fit parameters are presented in Table |B6| 


B2 XTE J0929-314: 

For ObsID 70096-03-02-00 (HF) data were extracted from 
PCUs 0, 2, 3 & 4 whereas for 70096-03-14-00 (LF) they 
were extracted from PCUs 0, 2, 4. The spectra were fit 
with a power-law continuum and thermal blackbody com¬ 
ponent. The hydrogen column density was fixed to Nh 


7.6 X 10^° cm ^ (Juett et al. 20031. The results are presented 


in Table 


B3 XTE J1807.4-294: 


For ObsIDs 70134-09-02-01 (HF) and 80145-01-17-02 (LF) 
the data were extracted from PCU2 whereas for 80419-01- 
01-01 (BE) data were extracted from PCUs 0, 2 and 3. The 


B7 HETE J1900.1-2455: 

The data were extracted from all layers of PCU2 for the 
considered ObsIDs. The h ydrogen column dens ity was fixed 
to Nh = 1 X 10^^ cm“^ (Papitto et al. 2013aI. We obtain 


the HF spectrum from ObsID 91015-01-06-00 with fit with 
tbabs* cutoff pi. We found a large excess at ~ 6 keV, which 
was modelled with a Gaussian centred at ~ 6.2 keV with 
a width of ~ 1 keV. The resulting fits were statistically 
acceptable. 

The ObsID 91059-03-02-00 (LF) was fit with 
tbabs (bbody+powerlaw). A weak Gaussian feature at 
~ 6.5 keV was added to improve the hts. The results are 
presented in Table [B7| 

A background estimate is not available for this source. 
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Table Bl. Spectral fit parameters of Swift J1756.9-2508. 



High flux 

Swift J1756.9-2508 
Low flux 

Background 

ObsIDs 

94065-02-01-05 

94065-06-02-03 

94065-06-03-02 

AfH(1022 cm-2) 

5.4 (fixed) 

5.4 (fixed) 

5.4 (fixed) 

Tbb (keV) 

- 

0.49 ± 0.08 

- 

Normbb 

- 

loti X 10-"! 

- 

r 

1.970 ±0.014 

2.24 ± 0.06 

2.67 ± 0.11 

X^/d.o.f. 

45.66/37 

40.39/41 

34.68/36 

Flux ( erg cm~^ 

6.30 ±0.04 X 10-1° 

1.990 ± 0.017 X 10-1° 

4.07 ± 0.15 X 10-11 


is the column density of neutral hydrogen for the tbabs model, Tr,b is the black body temperature, Norm^b is the normalisation of 
the bbody model, F is the power-law index. The first column represents the highest flux (HF) with pulsations, the second column gives 
the lowest flux (LF) with pulsations and the third column gives the observation used to measure the background emission (BE). 

Table B2. Spectral fit parameters of XTE J0929-314, see Table [BT] for details. 


XTE J0929-314 

High flux Low flux Background 


ObsIDs 

70096-03-02-00 

70096-03-14-00 

Ajy( 10°2 cm-2) 

0.0076 (fixed) 

0.0076 (fixed) 

Tbb (keV) 

0.82 ± 0.04 

0.53 ±0.05 

Normbb 

3.5 ±0.5 X 10-3 

1.9 ±0.4x10-° 

F 

1.80 ± 0.04 

1.90 ±0.11 

X^/d.o.f. 

36.33/42 

26.11/42 

Flux ( erg cm“^ s“^) 

4.42 ±0.02 X 10-1° 

6.64 ± 0.14 X 10-11 - 


B8 SAX J1808.4-3658: 

For ObsID 70080-01-01-04 (HF) the data were extracted 
from the top layers of PCU2. The hydrogen column den- 


sity was fixed to Nh = 2 x 10^^ cm ^ 

Papitto et al.||2009 

Cackett et al.||2009 Patruno et al. 2009c 1. The continuum 

was modelled with the thermal comptonisation NTHCOMP 

of 

Zdziarski et al. 

119961 

Zycki et al. ( 

19991. A blackbody 


and diskblackbody component were also required to model 
the continuum. A broad excess at ~ 6.2 keV was seen in the 
residuals, reminiscent of a relativistically broadened iron line 


reported in other works ( 

Papitto et al.||2009 Patruno et al. 

2009c Cackett et al. 2009 

I. However, owing to the poor spec- 


tral resolution of RXTE, we have not employed the sophis¬ 
ticated relativistic models like diskline and instead modelled 
the feature with a broadened Gaussian. The central energy 
of the Gaussian was fixed to 6.2 keV. 

Data for ObsIDs 30411-01-11-00 (LF) and 30411-01-11- 
02 (BE) was extracted from all PCUs and combined for 
spectral analysis. The data fit well with a simple power-law 
model. The results are presented in Table [BS] 


ponent. A narrow Gaussian component centred at ~ 6.5 keV 
(most likely a feature from background emission) was re¬ 
quired to obtain a statistically acceptable fit. The hydrogen 


column density was fixed to Nh = 2.87 x 10^^ cm ^ (|Torres 


et al.|20TT I. The results are presented in Table B9 


We note that flux measured from the LE and BE obser¬ 
vations are the same within error, indicating that the lowest 
flux with pulsations is background dominated. 


BIO XTE J1751-305: 


For all ObsIDs viz. 70131-03-01-00 (HF), 70131-01-09-000 
(LF) and 70131-02-04-00 (BE), data was extracted from 
PCU2. For the HF and BE, the spectra were fit with 
tbabs(po'werlaw+bbody). The LF spectrum was fit with a 
simple absorbed powerlaw. Additionally, a narrow Gaussian 
component centred at ~ 6.5 keV was required for the low 
flux observations. The neutral hydrogen column density was 
fixed to Nh = 1 x 10^^ cm“^ (Miller et al. 112003 Gierlifiski 


fc Poutanen|2005 1. The results are presented in Table BIO 


B9 IGR J1T498-2921: 

Data for all ObsIDs viz. 96435-01-02-01 (HF), 96435-01-06- 
04 (LF) and 96435-01-07-01 (BE) were extracted from all 
layers of PGU 2. All three spectra were fit with a simple 
power-law continuum model and a thermal blackbody com- 


Bll SAX J1748.9-2021: 

For ObsID 94315-01-06-07 (HF) and 60035-02-03-02 (LF) 
the data were extracted from PGU 2, whereas for 94315-01- 
11-02 (BE) the combined spectra from PGU 0 and 2 were 
used. The HP and LF spectra fit with tbabs(bbody+cutoffpl). 
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Table B3. Spectral fit parameters of XTE J1807.4—294, see Table [bT] for details. 


XTE J1807.4-294 

High flux Low flux Background 


ObsIDs 

70134-09-02-01 

80145-01-17-02 

80419-01-01-01 

Nh{10^'^ cm-2) 

0.56 (fixed) 

0.56 (fixed) 

0.56 (fixed) 

r 

1.90 ±0.03 

2.19 ±0.04 

2.16 ±0.04 

Tbb (keV) 

-‘-•^’^-0.19 

- 

- 

Normbb 

4.9 ±0.2 X 10“"^ 

- 

- 

X^/d.o.f. 

49.87/42 

31.07/44 

34.69/44 

Flux ( erg cm“^ 

8.19 ± 0.04 X 10-1° 

3.51 ± 0.07 X 10-1° 

7.25 ±0.15 X 10-11 


Table B4. Spectral fit parameters of NGC 6440 X-2, see Table [ bT] for details. 


NGC 6440 X-2 

High flux Low flux Background 


ObsIDs 

94044-04-02-00 

96326-01-40-01 

96326-01-36-00 

iVnClO^^ cm-2) 

0.59 (fixed) 

0.59 (fixed) 

0.59 (fixed) 

r 

1.83 ±0.02 

^■oD_o 12 

2.3 ±0.4 

r (Cutoflfpi) 

- 

1.1 ± 0.8 

- 

Ecut 

- 

5.5l° 

- 

X^/d.o.f. 

34.38/38 

45/35 

36.77/44 

Flux ( erg cm“^ s“^) 

2.62 ±0.03 X 10-1° 

3.4 ±0.2 X 10-11 

1.3 ±0.2 X 10-11 


Ecut is the high energy cutoff for a cutoff power-law model. 

Table B5. Spectral fit parameters of IGR J17511-3057, see Table [ bT] for details. 



High flux 

IGR J17511-3057 

Low flux 

Background 

ObsIDs 

94041-01-01-02 

94042-01-03-04 

94042-01-02-05 

Njj X 10^^ cm-^ 

1 (fixed) 

1 (fixed) 

1 (fixed) 

Tbb (keV) 

1.10± 0.16 

- 

- 

Nornibb 

6±3 X 10-1 

- 

- 

r 

1.70 ± 0.04 

2.02 ± 0.04 

2.37 ± 0.06 

X^/d.o.f. 

44.81/39 

39.89/41 

36.71/44 

Flux ( erg cm“^ s“^) 

8.65 ±0.05 X 10-1° 

1.000 ± 0.018 X 10-1° 

6.960 ± 0.018 X 10-11 


Table B6. Spectral fit parameters of XTE J1814-338, see Table [BT] for details. 



High flux 

XTE J1814-338 
Low flux 

Background 

ObsIDs 

80418-01-03-08 

80418-01-07-08 

80418-01-09-00 

Nh{10^'^ cm-2) 

0.167 (fixed) 

0.167 (fixed) 

0.167 (fixed) 

Tbb (keV) 

1.21 ± ±0.06 

- 

- 

Normbb 

1.7 ± 0.4 X 10-° 

- 

- 

r 

1.55 ± 0.04 

1.96 ±0.03 

2.4 ±0.3 

Ax^/d.o.f. 

63.99/42 

23.2/43 

28.17/44 

Flux ( erg cm“^ s“^) 

4.41 ± 0.03 X 10-1° 

6.00 ±0.10 X 10-11 

1.00 ±0.16 X 10-11 
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Table B7. Spectral fit parameters of HETE J1900.1—2455, see text and Table [ bT] for further details. 


HETE J1900.1-2455 

High flux Low flux Background 


ObsIDs 

91015-01-06-00 

91059-03-02-00 

cm-2) 

0.1 (fixed) 

0.1 (fixed) 

Tbb (keV) 


0.74 ± 0.10 

Nornibb 


8±2 X 10“"^ 

r (cutoffpi) 

0.9 ±0.1 


-Ecut 

4.0 ± 0.2 


r 


1.83 ±0.07 

X^/d.o.f. 

42.74/40 

25.26/40 

Flux ( erg cm“^ 

1.15 ±0.05 X 10-® 

3.84 ± 0.03 X 10 


Table B8. Spectral fit parameters of SAX J1808.4-3658. 



High flux 

SAX J1808.4-3658 
Low fiux 

Background 

ObsIDs 

70080-01-01-04 

30411-01-11-00 

30411-01-11-02 

cm-2) 

0.2 (fixed) 

0.2 (fixed) 

0.2 (fixed) 

r (Powerlaw) 

- 

2.24 ± 0.05 

2.23 ±0.11 

Tbb (keV) 

1.25 ± 0.07 

- 

- 

Nornibb 

9.2 ± 1.4 X 10-3 

- 

- 

Tdisk (keV) 

0.45 ± 0.06 

- 

- 

Nornidisk 

8.7)';5® X 103 

- 

- 

El (keV) 

6.2 (fixed) 

- 

- 

(Ti (keV) 

1.04 ±0.15 

- 

- 

r (NTHCOMP) 

2.3 ±0.1 

- 

- 

Te (keV) 

100 (fixed) 

- 

- 

X^/d.o.f. 

Flux ( erg cm“^ s“^) 

33 . 37/35 

1.850 ± 0.007 X 10“^ 

50.18/45 

2.82 ± 0.07 X 10-31 

35.47/43 

1.21 ± 0.08 X 10-11 


Tdisk is the temperature of the inner radius of the accretion disk, and Norm^isk is the normalisation of the diskbb model. The fluorescent 
iron line is fit with a Gaussian component, such that Ei is the line energy and <Ji the width of the line. Similarly E 2 and <J 2 are the 
line energy and width of a Gaussian component needed to model an excess at ~ 5.4 keV. See text and Table [BT] for further details. 

Table B9. Spectral fit parameters of IGR J17498—2921, see Table [bT] for details. 


IGR J17498-2921 

High flux Low flux Background 


ObsIDs 

96435-01-02-01 

96435-01-06-04 

96435-01-07-01 

cm-2) 

2.87 (fixed) 

2.87 (fixed) 

2.87 (fixed) 

Tbb (keV) 

1.51 ± 0.05 

1.93 ±0.15 

1.76 ±0.08 

Normbb 

1.9 ±0.2 X 10-3 

1.4ll°;3g X 10-3 

2.2 ±0.2 X 10-3 

r 

1.93 ±0.02 

2-39lg:j6 

2.3 ±0.1 

X^/d.o.f. 

52.76/39 

28.90/38 

32.39/38 

Flux ( erg cm~^ s~^) 

l'130l° °°3 X 10-9 

4.44 ± 0.04 X 10-19 

4.23 ±0.04 X 10-19 
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Table BIO. Spectral fit parameters of XTE J1751—305, see Table [BT] for details. 




XTE J1751-305 



High flux 

Low flux 

Background 

ObsIDs 

70134-03-01-00 

70131-01-09-000 

70131-02-04-00 

Aj^( 1022 cm-2) 

1 (fixed) 

1 (fixed) 

1 (fixed) 

Tbb (keV) 

1.95 ±0.16 

- 

- 

Nornibb 

2.1 ± 0.3 X 10-3 

- 

- 

Tdisk (keV) 

- 

- 

- 

Nornidiak 

- 

- 

- 

r 

1.77 ± 0.03 

2.070 ± 0.007 

2.32 ±0.03 

X^/d.o.f. 

47.59/40 

38.07/35 

56.08/40 

Flux ( erg cm“^ 

1.500to;o'(® ^ 10 "® 

3.97 ± 0.01 X 10-3° 

6.51 ± 0.08 X 10-13 


Table Bll. Spectral fit parameters of SAX J1748.9—2021, see Table [BT] for details. 




SAX J1748.9-2021 



High flux 

Low flux 

Background 

ObsIDs 

94315-01-06-07 

60035-02-03-02 

94315-01-11-02 

Ah( 1032 cm-2) 

0.59 (fixed) 

0.59 (fixed) 

0.59 (fixed) 

Tbb (keV) 

2.11 ± 0.14 

O-dSto.?! 

- 

Normbb 

1.30 ±0.19 X 10-3 

X 10-3 

- 

r (cutoffpi) 

1.00 ±0.16 

^•'jy_0.23 

- 

Ecut (keV) 

4.1 ± 0.6 


- 

r (power law) 

- 

- 

2.2 ± 0.2 

xVd.o.f. 

45.41/41 

36.21/33 

28.23/43 

Flux ( erg cm~^ s~^) 

4.13 ±0.01 X 10-® 

2.960 ±0.005 X 10-® 

1.9 ±0.2 X 10-11 


The BE spectrum was fit using an absorbed power-law. The 
hydrogen column density was fixed to Nh = 5.9 x 10^^ cm“^ 
| |Harris|[l996[ ). A Gaussian emission feature centred at ~ 
6.5 keV was used to improve statistics of the LF spectrum 
ht. The results are presented in Table [BTTI 


B12 Swift J1749.4-2807: 


The data were extracted from PCU2 for all ObsIDs (HF: 
95085-09-01-00, LF: 95085-09-02-07, BE: 95085-09-02-10). 
There was an eclipse during the observation in 95085-09-01- 


00 (Markwardt & Strohmayer 2010 Ferrigno et al. 20111, 


so for the spectral analysis of that observation we consid¬ 
ered only the data from the initial 990 s when the source 
was visible. The spectra were fit with power-law continuum 
model. The neutral h ydrogen column density was fixed to 
Nh = 3 X 10^^ cm“^ i Ferilgno et al.|r2011 Wijnands et al. 
2009). For the LF and BE an additional blackbody compo¬ 


nent was required. The results are presented in Table |BT^ 


B13 Aql X-1: 

Since the pulsations of Aql X-1 have been detected only once 
(ObsID 30188-03-05-00), we use this observation as both 
the HF and LF. We extract the data from PCU2 and ht 


the spectrum with a tbabs(cutoffpl+Gaussian+bbody) model. 
The background emission is evaluated from ObsID 30073-06- 
01-00, which was ht with a simple absorbed power-law. Our 
results are shown in Table IbT^ 


B14 IGR J00291+5934: 


For the ObsID 90052-03-01-00 (HF) data from PCU2 were 
used for spectral analysis. For 90425-01-02-01 (LF) and 
90425-01-03-06 (BE) data from all active PCUs (PCU 0, 
2 and PCU 0, 2, 3 respectively) were combined to improve 
the photon statistics. 

The HF spectrum was ht with the thermal blackbody 
and power-law models. The LF and BE spectra were ht 
with a simple absorbed power law. The neutral hydrogen 
column density was hxed to Nh ~ 0.46 x 10^^ cm“^ based 
on measurements from XMM-Newton and Chandra obser¬ 


vations (Torres et al. 2008 Paizis et al. 2005). A gaussian 


emission feature was required to improve the hts of the LF 
and BE spectra. The ht results are given in Table [Bid] 
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Table B12. Spectral fit parameters of Swift J1749.4-2807, see Table [ bT] for details. 



High flux 

Swift J1749.4-2807 

Low flux 

Background 

ObsIDs 

95085-09-01-00 

95085-09-02-05 

95085-09-02-10 

Nh{10^'^ cm-2) 

3 (fixed) 

3 (fixed) 

3 (fixed) 

Tbb (keV) 

- 

1.67 ± 0.11 

1.8 ±0.2 

Normbb 

- 

7.94 ± 1.2 X 10-1 

5.8 ± 1.2 X 10-1 

r 

1.89 ±0.02 

1.78 ±0.08 

1.90 ±0.09 

X^/d.o.f. 

43.21/41 

54.01/33 

39.77/39 

Flux ( erg cm“^ 

5.24 ±0.06 X 10-1° 

2.67 ± 0.02 X 10-1° 

2.41 ± 0.03 X 10-1° 


Table B13. Spectral fit parameters of Aql X-1, see Table [ bT] for details. 



High/Low flux 

Aql X-1 

Background 

ObsIDs 

30188-03-05-00 

30073-06-01-00 

NnilO'^^ cm-°) 

0.4 (fixed) 

0.4 (fixed) 

Tbb (keV) 

+0.44 

^•^-*--0.09 

- 

Normbb 

-‘-•^-0.6 

- 

r 

- 

1.92 ±0.02 

F (cutoffpl) 

6 6 

1 

r- 

o 

- 

Tcut (keV) 

3.4 ± 0.4 

- 

El (keV) 

o.az_o 38 

- 

CTi (keV) 

1.0 ±0.3 

- 

X°/d.o.f. 

45.95/38 

45.39/43 

Flux ( erg cm~^ s“^) 

8.74 ±0.01 X 10 

-° 1.34 ± 0.05 X 10-11 


Table B14. Fit parameters from spectral analysis of IGR J00291+5934. 



High flux 

IGR J00291+5934 
Low flux 

Background 

ObsIDs 

90052-03-01-00 

90425-01-02-01 

90425-01-03-06 

IV^f(1022 cm-2) 

0.46 (fixed) 

0.46 (fixed) 

0.46 (fixed) 

r 

1.50 ±0.04 

1.73 ±0.04 

1.50 ±0.03 

Tbb (keV) 

1.14 ± 0.06 

- 

- 

Normbb 

1.05 ±0.03 X 10-° 

- 

- 

X°/d.o.f. 

48.81/42 

22.92/35 

30.55/41 

Flux ( erg cm~^ s“^) 

9.7 ± 0.05 X 10-1° 

1.09 ±0.02 X 10-1° 

5.76 ±0.09 X 10-11 











